Magnetic molecules created by hydrogenation of Poly cyclic Aromatic Hydrocarbons 
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Present routes to produce magnetic organic-based materials adopt a common strategy: the use of magnetic 
species (atoms, polyradicals, etc.) as building blocks. We explore an alternative approach which consists of 
selective hydrogenation of Polycyclic Aromatic Hydrocarbons. Self-Consistent-Field (SCF) (Hartree-Fock and 
DFT) and multi-configurational (CISD and MCSCF) calculations on coronene and corannulene, both hexa- 
hydrogenated, show that the formation of stable high spin species is possible. The spin of the ground states 
is discussed in terms of the Hund rule and Lieb's theorem for bipartite lattices (alternant hydrocarbons in this 
case). This proposal opens a new door to magnetism in the organic world. 

PACS numbers: 31.10.+Z, 33.15.Kr, 31.15.aq, 75.50.Xx 



I. INTRODUCTION 

Two successful routes that are being actually followed to 
produce magnetic organic materials^ are the addition of mag- 
netic atoms^ and the use of polyradicals^. In particular, 
carbon-based nickel compounds that show spontaneous field- 
dependent magnetization and hysteresis at room temperature, 
have been recently synthesized^. Moreover, the combination 
of two radical modules with different spins has allowed the 
obtaining of organic polymers with ferro- or antiferromag- 
netic ordering"^. Research on molecules containing polyrad- 
icals goes back to the early nineties^ ^ ^ ^ and has produced a 
variety of results as, for example, the synthesis of high spin 
organic molecules. In some of these molecules the failure 
of Hund's rule has been demonstrated^. On the other hand, 
experimental and theoretical evidence has been recently pre- 
sented indicating that 5-dehydro-m-xylylene or DMX was the 
first example of an organic tri-radical with an open- shell dou- 
blet ground-statd^. Both methods share a common strategy: 
the use of ingredients (either radicals or atoms) that provide a 
finite spin. 

In this work we follow a different approach. Specifically, 
we predict the existence of spin polarized organic molecules 
derived from non magnetic 7r-conjugated Polycyclic Aromatic 
Hydrocarbons (PAHs) by selective hydrogenation of their pe- 
ripheral C atoms. High hydrogenation of PAHs has been pro- 
posed as a method for hydrogen storage^. More recently, the 
feasibility of double hydrogenation of those compounds has 
been investigated theoretically^. 

Our work is inspired upon Lieb's theorem for bipartite lat- 
tices that shows the appearance of magnetism whenever they 
are unbalanced^^. According to Lieb, if a nearest neigh- 
bor model with a local on-site interaction is applicable to a 
bipartite lattice, the spin multiplicity of the ground state is 



I Na — Nb 1 + 1, where Na and Nb are the number of atoms in 
each sublattice. Most PAHs are alternant hydrocarbons where 
carbon atoms can be separated into two disjoint subsets so that 
an atom in one set only has neighbors in the other set (Figs. 1 
and 2 show a colored version of the partition). The same the- 
orem has been used to support the existence of magnetism in 
graphene ribbons and islands^. All work we know is based on 
single-determinantal methods, i.e., on a more or less sophis- 
ticated form of Self-Consistent-Field (SCF) calculation. Let 
us remark that being 7r-orbital magnetism a direct result of the 
strong correlation among 7r-electrons, only methods designed 
explicitly to catch these effects (like CISD and MCSCF, used 
in our work) can help to resolve the doubts regarding the ap- 
pearance of magnetism in graphite-derived systems. 

The rest of the paper is organized as follows. The ab ini- 
tio methods (both mono- and multi-determinantal) used in this 
work are discussed in some detail in section II, while the re- 
sults obtained with those methods are reported and discussed 
in section III. Section IV in turn is devoted to the analysis of 
the ab initio results by means of model Hamiltonians, in par- 
ticular the Hubbard and the Pariser-Parr-Pople Hamiltonians. 
Finally, the conclusions of our work are gathered in section V. 



II. AB INITIO CALCULATIONS: METHODS AND 
NUMERICAL PROCEDURES 

Calculations of the spin states of the molecules of Figs. 
1 and 2 were done using the following basis functions sets: 
MIdP, cc-pVDZ and cc-pVTzP. Although the latter set 
guarantees a sufficient precision, varying the dimension of 
the variational space allowed to check the reliability of our 
results. SCF calculations were carried out at the Restricted- 
Hartree-Fock (RHF) level and by means of the hybrid den- 





FIG. 1: (color online) 1, 4,5, 8,9, 12-hexahydrocoronene (A, hereafter 
referred to as Ds^ according to its symmetry group), 1,3,5,7,9,11- 
hexahydrocoronene (B, hereafter referred to as Cs/i) and planar 
1,4,5,6,7, 10-hexahydrocorannulene (C, hereafter referred to as C2v)- 
Saturated carbon atoms are represented by black symbols while dark 
gray (magenta) and light gray symbols are used to distinguish car- 
bon atoms belonging to different sublattices. Corannulene is a non- 
alternant hydrocarbon, that is, a frustrated cluster of carbon atoms 
(note the fully magenta bond between two magenta atoms). 




TABLE I: Total energies (in Hartrees) for atomic hydrogen, molecu- 
lar hydrogen, coronene C24H12, corannulene C20H10, two molecules 
obtained from hexahydrogenation of coronene and one derived from 
hexahydrogenation of corannulene (in the latter case the results cor- 
respond to the planar geometry shown in Fig. IC). The results were 
obtained using three basis sets (MIDI, cc-pVDZ and cc-pVTZ), two 
SCF methods (RHF and RB3LYP) and one multi-configurational 
method CISD. The number of occupied (m) and empty (n) tv Molecu- 
lar Orbitals included in the CISD calculations as well as the number 
of electrons that fill them (N) is indicated as (m-i-n,N). Small stars 
emphasize the spin multiplicity of the more stable state. 



MOLECULE 


METHOD 


BASIS SET 


MIDI cc-pVDZ cc-pVTZ 


H 


RHF 


-0.4970 -0.4993 -0.4998 




RB3LYP 


-0.4953 -0.4979 -0.4988 


H2 


RHF 


-1.1217 -1.1287 -1.1330 




RB3LYP 


-1.1623 -1.1668 -1.1733 


Coronene 


RHF 


-910.4869 -916.0197 -916.2293 


C24H12 


RB3LYP 


-915.9341 -921.3874 -921.6253 


Corannulene 


RHF 


-758.6127 -763.2326 -763.4078 


C20H10 


RB3LYP 


-763.1633 -767.7138 -767.9092 


C24H18 


RHF 


-913.6339 -919.2239 -919.4337 


^ D3h (S=3) ^ 


RB3LYP 


-919.1814 -924.6701 -924.9159 




CISD(1 1-^3,16) 


-913.7112 -919.2905 -919.4939 


C24H18 


RHF 


-913.3714 -918.9826 -919.2040 


Dsh (S=0) 


RB3LYP 


-919.0751 -924.5639 -924.8175 




CISD(8-h6,16) 


-913.4988 -919.0286 -919.2431 


C24H18 


RHF 


-913.5690 -919.1638 -919.3742 


Csh (S=3) 


RB3LYP 


-919.1228 -924.6146 -924.8620 




CISD(1 1-^3,16) 


-913.6153 -919.2027 -919.4078 


C24H18 


RHF 


-913.7732 -919.3607 -919.5746 


^ Csh (S=0) ^ 


RB3LYP 


-919.3423 -924.8230 -925.0730 




CISD(8-h6,16) 


-913.8433 -919.4158 -919.6110 


C20H16 


RHF 


-761.9008 -766.5625 -766.7408 


• C2v (S=2) • 


RB3LYP 


-766.5503 -771.1250 -771.3348 




CISD(8-h4,12) 


-761.9652 -766.6291 -766.8014 


C20H16 


RHF 


-761.7591 -766.4418 -766.6266 


C2v (S=0) 


RB3LYP 


-766.4913 -771.0760 -771.2842 




CISD(6-h6,12) 


-761.8264 -766.4643 -766.6588 



FIG. 2: (color online) Two views of curved 1,4,5, 6,7, 10- 
hexahydrocorannulene in the calculated stable geometry (hereafter 
referred to as Ci). As in Fig. 1, black symbols indicate carbon atoms 
forming only single bonds while dark gray (magenta) and light gray 
symbols denote each of the two sublattices in which carbon atoms 
can be separated. 



sity functional RBSLYfIHEIIII. in both cases the Restricted- 
Open- Shell variant was used in order to get well-defined to- 
tal spin values ^^. In order to check the accuracy of the de- 
scription of the correlation energy of partially filled tt- shells, 
multi-configurational wave-functions calculations were also 
performed. Configuration Interaction with Single and Dou- 



ble excitations (CISD) calculations^ were carried out in all 
cases, while some checks were also made by means of the 
Multi-Configurational SCF (MCSCF) on the fully optimized 
set in the active space version^^^^. The active space was gen- 
erated within the following windows (m-Fn,N) of m occupied 
and n empty tt Molecular Orbitals (MO) filled with N elec- 
trons: hexahydrogenated coronene S=0, (8+6,16) and S=3, 
(11+3,16), and planar hexahydrogenated corannulene S=0, 
(6+6,12) and S=2, (8+4,12). Other tt-MO lie excessively far 
from the HOMO-LUMO gap to give a sizable contribution. 
Geometries were only optimized at the SCF (RB3LYP) level. 
The geometry of 6H-corannulene was optimized for both its 
planar metastable form and its curved stable form (see Figs. 




TABLE II: Fragmentation energies (in Hartrees) of molecules de- 
rived from hexahydrogenation of coronene and of corannulene. To- 
tal energy differences are given both for atomic and molecular forms 
of hydrogen. Data of Table I have been used and again small stars 
emphasize the spin multiplicity of the more stable state. 



FIG. 3: (color online) Total spin densities of 1,4,5,8,9,12- 
hexahydrocoronene and 1,3,5,7,9,1 1-hexahydrocoronene (both 
corresponding to septuplets, S=3) and planar 1,4,5,6,7,10- 
hexahydrocorannulene (S=2 state). 



IC and 2). However, in order to allow a discussion in terms of 
TT-orbital models, the results for the energies of its spin states 
discussed hereafter correspond to the planar geometry. Any- 
how, energy differences between the spin states of the two 
allotropes are very small (fragmentation energies for both pla- 
nar and curved geometries are reported below). All quan- 
tum chemistry calculations were done using the GAMES S 
progranPi]. 



III. AB INITIO CALCULATIONS: RESULTS 

Total energies for the singlet and the relevant multiplet of 
hydrogenated coronene Dsh, Csh and planar hydrogenated 
corannulene C2v (A, B and C in Fig. 1) are reported in Ta- 
ble I. It is first noted that whereas the energies obtained with 







Atomic H 


Molecular H2 


MOLECULE 


METHOD 


MIDI cc-pVTZ 


MIDI cc-pVTZ 


C24H18 

• i:>3^(s=3)^ 


RHF 
RB3LYP 


-0.1651 -0.2056 
-0.2755 -0.2980 


0.2181 0.1946 
0.2396 0.2292 


C24H18 

Dsh (S=0) 


RHF 
RB3LYP 


0.0974 0.0242 
-0.1692 -0.1996 


0.4805 0.4243 
0.3459 0.3276 


C24H18 

C3h (S=3) 


RHF 
RB3LYP 


-0.1002 -0.1460 
-0.2170 -0.2441 


0.2830 0.2542 
0.2982 0.2831 


C24H18 

• C3h (S=0) ^ 


RHF 
RB3LYP 


-0.3044 -0.3465 
-0.4365 -0.4551 


0.0788 0.0537 
0.0787 0.0720 


C20H16 
• C2v (S=2) ^ 


RHF 
RB3LYP 


-0.3061 -0.3342 
-0.4152 -0.4290 


0.0770 0.0660 
0.0999 0.0981 


C20H16 
C2v (S=0) 


RHF 
RB3LYP 


_0.1644 -0.2200 
-0.3562 -0.3784 


0.2187 0.1802 
0.1589 0.1487 


C20H16 
ic Ci (8=2) • 


RHF 
RB3LYP 


-0.3024 -0.3259 
-0.4098 -0.4198 


0.0807 0.0743 
0.1054 0.1073 


C20H16 
Ci (S=0) 


RHF 
RB3LYP 


-0.1607 -0.2072 
-0.3471 -0.3643 


0.2224 0.1930 
0.1680 0.1628 



the small basis set MIDI and those obtained with the already 
large cc-pVDZ, differ in 4-6 Hartrees (approximately 0.6%), 
the difference is reduced to 0.1-0.3 Hartrees (approximately 
0.02%) when cc-pVDZ is replaced by the largest basis used 
in this work, namely, the cc-pVTZ basis set. This indicates 
that convergence, as far as the basis set is concerned, is rather 
acceptable. In the case of hexahydrogenated coronene (briefly 
6H-coronene), results clearly show that, no matter the method 
or the basis set used, the ground state of molecule Dsh is a 
septuplet and that of molecule Csh a singlet. We have checked 
that other spin states lie between those two. In molecule Dsh 
the largest energy difference between the high spin ground 
state and the singlet occurs for RHF (0.23-0.26 Hartrees). 
This difference is reduced to approximately 0.1 Hartrees for 
RB3LYP, increasing again using the CISD method. On the 
other hand, all results for Csh conformation show that the sin- 
glet is below the septuplet by more than 0.2 Hartrees. Simi- 
lar results are obtained for 6H-corannulene, although energy 
differences are slightly smaller. Table I also reports total en- 
ergy results for atomic and molecular hydrogen, coronene and 
corannulene that allow the calculation of fragmentation ener- 
gies (Table II analysis). These are negative relative to atomic 
hydrogen but not relative to the molecular form. Therefore, 
actual synthesis of the hydrogenated molecules would need 
sophisticated reaction paths^^. We also note that the sin- 
glet ground state of Csh hydrogenated coronene is more sta- 
ble than that of the molecule having a septuplet ground state 
(Dsh)' Presumably, other forms of 6H-corannulene would 
also show deeper ground state energies than that of the stud- 
ied magnetic conformation. Note also that hydrogenation of 



the curved (stable) geometry of corannulene (see Fig. 2) is 
slightly less favorable than that of its planar geometry (com- 
pare results for C2v and Ci in Table II). Anyhow, as in the 
planar geometry, the quintuplet has a lower energy than the 
singlet. 

Fig. 3 depicts the total spin densities of the septuplet states 
(S=3) of 1,4,5,8,9,12-hexahydrocoronene and 1,3,5,7,9,11- 
hexahydrocoronene (A and B) and the quintuplet (S=2) 
of planar 1,4,5,6,7,10-hexahydrocorannulene. Concerning 
1,4,5,8,9,12-hexahydrocoronene, the most appealing result is 
that the spin density is finite only on the carbon atoms of one 
sublattice. More precisely, spin density is located in the sub- 
lattice to which no additional H atoms were attached. This 
result is highly illustrative allowing some intuition on the rea- 
sons for a magnetic ground state: electron-electron repulsion 
is minimized because each electron avoids sitting at nearest- 
neighbors distances from the others. However, in 1,3,5,7,9,1 1- 
hexahydrocoronene, a molecule with a singlet ground state, 
the spin is equally spread over the two sublattices implying 
larger electronic repulsions at the central hexagon. The case 
of 1,4,5,6,7, 10-hexahydrocorannulene is even more interest- 
ing as, being a frustrated molecule, at least one bond be- 
tween atoms of the same sublattice should be present. This 
is clearly visible in Fig. 3 once a sublattice is identified as 
the sites showing spin density while the rest belong to the 
other sublattice (Colors in Fig. 1 have anticipated this fea- 
ture). We will show later that the model Hamiltonian calcula- 
tions for 1,4, 5, 6,7, 10-hexahydrocorannulene show frustration 
at the same bond than ah initio calculations (compare Figs. 3 
and 4). Having identified the atoms at each sublattice, it is 
tempting to use the unbalance in the molecule {Na — Nb=^) 
to predict the total spin of the ground state using Lieb's for- 
mula. The result (S=2) is in perfect agreement with numerical 
results. This is particularly interesting as in principle Lieb's 
theorem should only work on non-frustrated systems. 

Spin multiplicity of the ground state of a molecule is usu- 
ally predicted by means of Hund rule applied to MO energies 
obtained by an appropriate method. We have checked that the 
spin of the ground states of the molecules here investigated 
is consistent with the degeneracy of the HOMO that Hiickel's 
method gives for the skeleton of C atoms having an unsatu- 
rated TT orbital. This is true not only for 6H-coronene, but also 
for 6H-corannulene. Although the extended Hiickel's method 
used by ab initio codes to initialize the self-consistency pro- 
cess slightly lifts this degeneracy, the HOMO still appears as a 
narrow bunch containing a number of orbitals compatible with 
the spin of the ground states of the three planar molecules de- 
picted in Fig. 1. Then, as in Hund rule, such a distribution 
of molecular orbitals favors high spin ground states through 
a winning competition of interaction energy gains against ki- 
netic energy losses. 



IV. MODEL HAMILTONIANS 

Let us critically examine the applicability of Lieb's theo- 
rem as the predicting tool of the multiplicity of the ground 
state of hydrogenated PAHs. The underlying Hubbard model 



ignores that: (i) transfer integrals in any realistic system are 
not limited to nearest neighbors sites, (ii) cr-orbitals appear 
around the HOMO-LUMO gap in the same energy interval 
as TT-orbitals, (iii) interaction among electrons is not limited 
to on-site Coulomb repulsion. In our opinion, the success of a 
theorem or rule based on the simplest interacting model comes 
from its actual capability of describing the correct antiferro- 
magnetic spin-spin correlations between nearest n electrons. 
Strong correlation is the basis for the basic correctness of a 
simplified image in which up and down spins altemate^^ . 

Even if the spin multiplicity of the ground state is predicted 
either by Hund rule or Lieb's theorem, a deeper understand- 
ing of underlying correlations calls for a complete numerical 
solution of simple interacting models. We have analyzed both 
Pariser-Parr-Pople (PPP) model Hamiltonian^"^ ^^ and the local 
version of Hubbard Hamiltonian^^, which actually is a partic- 
ular case of the former. The PPP Hamiltonian contains a non- 
interacting part ^0 and a term that incorporates the electron- 
electron interactions Hj : 



The non-interacting term is written as. 



Ho 



^0 E 

i=l,N;a 



C-ri^Cj. 



t 



<ij>\(7 



(1) 



(2) 



where the operator c\^ creates an electron at site i with spin a, 
eo is the energy of carbon 7r-orbital, and t is the hopping be- 
tween nearest neighbor pairs (kinetic energy). N is the num- 
ber of unsaturated C atoms. The interacting part is in turn 
given by. 



Hi = U ^ rii^riii + ^ E ^i-jl'^i^jrija' (3) 



i=l,N;a 



iy^j;a,a' 



where U is the on-site Coulomb repulsion and V\i_j\ is the 
inter-site Coulomb repulsion, while the density operator is 



-.t 



^icr^icr • 



(4) 



This Hamiltonian reduces to the Hubbard model for Vj^_j|=0 
We start discussing the fitting of spin state energies by 
means of Hubbard Hamiltonian. Lanczos algorithm in the 
whole Hilbert occupation space is used to get numerically 
exact many -body states^^. Coulomb on-site repulsion has 
been adjusted to describe spin excitations of some PAHs. Be- 
cause the interacting model cannot be solved exactly for 6H- 
coronene (18 orbitals or equivalently sites lead to a Hilbert oc- 
cupation space of dimension equal to 4^^ which is beyond ac- 
tual computational facilities), the case of a coronene molecule 
with all peripheral C atoms saturated by additional H has been 
considered. This leaves a molecule with only 12 tt orbitals, a 
cluster size that can be easily handled by means of Lanczos al- 
gorithm. Also benzene (6 sites), anthracene (14 sites) and 6H- 
corannulene (14 sites) have been fitted. Calculations were car- 
ried out by taking the hopping integral commonly used to de- 
scribe graphene sheets, t = — 2.71 eV, and varying the on-site 
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FIG. 4: (color online) Left: Total RB3LYP (circles) and Hubbard model (continuous curves) energies of spin S states versus the on-site 
Coulomb repulsion U for dodecahydrogenated coronene C24H24 (each peripheral carbon atom saturated with an additional hydrogen), an- 
thracene and 6H-corannulene C20H16 (Fig. IC). The sequence from lowest to highest energy is S=0,1, ... for anthracene and C24H24 and 
8=2,1,0,3,4, ... for C20H16 (in this molecule the Hubbard model gives an almost twofold degenerate ground state). Energies are referred to 
the corresponding state of maximum spin, except for C24H24 that were downward shifted by 1.8 eV to improve the overall fit (the state of S=6 
is largely participated by H orbitals turning invalid the Hubbard model). Molecular geometries were only optimized for the ground state and 
taken unchanged for the calculation of excited spin states. Transfer integral was taken equal to -2.71 eV, and, as the results indicate, U = 3.3 eV 
nicely reproduces the ab initio energies. Right up: Spin-spin correlations in C20H16 calculated by means of the Hubbard model (the skeleton 
of C atoms is shown as an inset) on blue pair (two atoms placed on the symmetry axis) and red pair (frustrated horizontal bond). Right down: 
Ab initio ground state energies (circles) of the charged states of a benzene molecule that is not allowed to relax. Energy differences are plotted 
relative to the neutral case. Results obtained by means of Hubbard (triangles) and Pariser-Parr-Pople (squares) models are also shown. 



repulsion U. The results depicted in the left panel of Fig. 4 in- 
dicate that spin states of these molecules can be reasonably fit- 
ted with /7 = 3.3 eV (benzene, for which the fitting is as good 
as for anthracene, is not shown for the sake of clarity). Albeit 
noticeable deviations occur in the three lowest lying states of 
6H-corannulene, the state ordering is the correct one. We have 
checked that the failure to correctly separate the lower excita- 
tions of 6H-corannulene is not exclusive of the simplest in- 
teracting model: a PPP calculation using Ohno's interpolation 
scheme^^ shows a similar weakness. Let us remark once more 
that, despite the rather small U (U/\t\ = 1.27) resulting from 
the fittings shown in Fig. 4, anti-ferromagnetic correlations 
in these molecules, as calculated by means of the interacting 
model, are significant. Particularly attracting is the case of 
corannulene for which there is a bond at which the spin- spin 
correlation is significantly smaller than at other bonds of the 
molecule (top right panel of Fig. 4). Interestingly enough, 
placing frustration (two adjacent tt orbitals of the same sub- 
lattice as shown in Fig. IC) at that bond gives a difference 
between carbon atoms in the two sublattices of four atoms, 
which, using Lieb's formula, predicts a ground state of total 
spin 2, in agreement with our numerical results. Summariz- 
ing, U = 33 eV works satisfactorily describing the spin states 
of aromatic molecules and, in particular, the multiplicity of 
the ground state. 



The same simple Hubbard model fails, however, in describ- 
ing the charged states of these systems as results for benzene 
show (bottom right panel of Fig. 2). Lanczos results for the in- 
teracting Hamiltonian are compared with B3LYP results for a 
charged benzene ideally restricted to a fixed geometric struc- 
ture. Actual energy differences are much higher than those 
predicted by the model. Our results do also illustrate the lack 
of electron-hole symmetry that characterizes any realistic self- 
consistent field calculation as opposed to Lieb's model. The 
PPP model with t = —2.71 eV and values for the Coulomb 
repulsion integrals from Ref. |^ although greatly improves 
the fitting, still gives a symmetric curve by the well-known 
pairing between occupied and unoccupied molecular orbitals. 
A full fitting of charge and spin states will surely require in- 
cluding a larger number of parameter d^^'^^l 



V. CONCLUDING REMARKS 

We have proposed a new route to produce magnetic or- 
ganic molecules that consists of hydrogenating PAHs. In the 
case of alternant PAHs the spin multiplicities of the molecules 
ground states agree with Lieb's prediction, even though ab 
initio Hamiltonians may significantly differ from Hubbard's 
model. A probably related result is that ab initio energies of 



the spin states of these molecules can be very satisfactorily 
fitted by means of the simplest version of the Hubbard model. 
It seems that the molecule topology is enough support for the 
main result. Energies of charged molecules, instead, cannot 
be described by the simple, most popular, interacting models, 
suggesting a critical examination of their use in, for instance, 
graphene. Results for total spin densities in molecules having 
a magnetic ground state clearly show that the spin is localized 
in only one of the two sublattices. On the other hand, ab initio 
and model Hamiltonian calculations for hydrogenated coran- 
nulene, place the frustrated bond at the same location. This 
produces an unbalance in the molecule which, using Lieb's 
formula, gives S=2 for the ground state of the molecule, in 
agreement with the numerical results. It is also worth noting 
that in a recent stud>^^ we have shown that dehydrogenation 
may also produce magnetic molecules. Although dehydro- 
genation is a highly unlikely process, dehydrogenated PAHs 



have been intensively investigated by astrophysicists^^ who 
believe they to form part of interstellar matter. Although the 
results presented here are encouraging, there is still a long way 
to go: finding procedures to synthesize these hydrogenated 
PAHs and crystallize them into solids that may eventually 
show magnetic properties. 
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